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We consider several effects of the matter wave dynamics which can be observed in Bose- 
Einstein condensates embedded into optical lattices. For low-density condensates we 
derive approximate evolution equations, the form of which depends on relation among 
the main spatial scales of the system. Reduction of the Gross-Pitaevskii equation to a 
lattice model (the tight-binding approximation) is also presented. Within the framework 
of the obtained models we consider modulational instability of the condensate, solitary 
and periodic matter waves, paying special attention to different limits of the solutions, 
i.e. to smooth movable gap solitons and to strongly localized discrete modes. We also 
discuss how the Feshbach resonance, a linear force, and lattice defects affect the nonlinear 
matter waves. 



1. Introduction 

The first experimental realization of the Bose-Einstein condensate"'^ (BEC) resulted 
in explosion of theoretical and experimental research devoted to the new state of 
matter, which now attracts a great deal of attention of the physical community. 
The studies have become really interdisciplinary, involving scientists working not 
only in atomic physics, but also in nonlinear optics, condensed matter physics, 
dynamical systems, differential equations, etc. Analogies between the theory of 
nonlinear matter waves in BECs and earlier developments in other areas of the 
nonlinear science significantly accelerate the research progress. The major part of 
the knowledge available today is gathered in a number of recent reviews^. 

It is widely accepted that generation, dynamics and management of nonlinear 
matter waves are the most interesting and important issues of the mean-field the- 
ory of a BEC. In particular, three types of relatively stable coherent structures 
play the main role in the theory. They are bright solitons, characterized by atomic 
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density distribution localized in space; dark solitons, representing holes in the oth- 
erwise homogeneous atomic distributions; and nonlinear periodic waves. Although 

these excitations can be observed experimentally in BECs trapped by a magnetic 
field^'^'^, due to the experiment^ it became clear that the matter wave dynamics 
can be drastically changed and appreciably enriched if an optical lattice (OL) is 
imposed on the condensate (in addition to or instead of a harmonic trap). That 
is why BECs loaded into OLs, which are created by laser beams, attract especial 
attention. A large diversity of remarkable effects has already been observed ex- 
perimentally. Among them we mention Bloch oscillations of BECs®'^, instability 
of nonlinear matter waves^, Landau-Zcncr tunneling^, supcrfluid to Mott-insulator 
phase transition^", compression of a condensate^ ^, and gap matter solitons^^. 

Behavior of a BEC strongly depends on parameters of an OL which in its turn 
can be controlled by intensities and/or by geometry of laser beams. In other words, 
the matter wave dynamics can bo effectively manipulated by external fields, the 
problem having primary importance for a number of BEC applications, like atomic 
interferometry^^, quantum phase transitions^^, matter wave amplification^^, atom 
laser^^ and others. 

The scope of this brief review does not allow us to cover all important issues. 
Instead, we deal with a few topics of the mean-field theory of BECs embedded in 
OLs, only mentioning, in Conclusion, other important issues. 



2. Mean-filed approximation 

Quantum description of a BEC in the approximation of two-body interactions is 
provided by the nonlinear Schrodinger (NLS) equation for the field operators = 
*(r,t)i7: 



ih- 



dt 



where Vext{^) is an external trap potential, A4 is an atomic mass, and ^(r) de- 
scribes interaction between two atoms. In the present paper the consideration will 
be restricted to local interactions, when one can approximate V(r) — 5o^(r) with 
go = ATrfi^as/M. Moreover, at zero temperature, the operator ^ can be repre- 
sented in the form ^ = + ip, where is an algebraic function, called an order 
parameter (or macroscopic wave function) and the operator ip describes quantum 
fluctuations^^. Then, neglecting quantum fluctuations^^, the dynamics of the con- 
densate can be described by the macroscopic wave function which solves the 
Gross-Pitaevskii (CP) equation^^ 

ih— = -—A^ + Ve.t{r)^ + go\^f^ (2) 
and is normalized to the number of condensed atoms J\f: 

\'^fd^r=Af. (3) 



I 



I 
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In the case of a BEC embedded in an OL the external potential usually consists 
of two parts: Vext{^) — Vt{j^) + ^^(r). The first term is a harmonic trap given by 
Vt(i") = ^ + ti^o^^) where r = (a;,r_L), = {y,z), and uij_ and are a 

transverse and axial harmonic oscillator frequencies (for the sake of simplicity the 
consideration will be restricted to radially symmetric traps). Respectively, one can 
introduce two spatial scales a± = y^h/'moj± and qq = y^h/mujQ, characterizing 
transverse and longitudinal extensions of the atomic cloud. The second term is a 
periodic potential which, almost everywhere through this paper is considered to be 
one-dimensional (ID), VL{r) = Vl{x) = Vl{x + A), with a period A. 

Interaction of a single atom, having an electric dipole moment d, with an electric 
field E is described by the Hamiltonian Hd = — d • E, which results in the shift of 
the energy level (see e.g.^°): ^^(r) — — ^a'(E^)t where (•)( stands for time averaging 
and a' is the real part of the polarizability. Let us consider two monochromatic laser 
beams, E — eii?^e"^^* + e2i5_c^e~*"^*, where ujl is a frequency, are polarization 
verctors, and — E^^ (an overbar stands for the complex conjugation), in a 
geometry^^ typical for many experimental settings^ and shown in Fig. 1, Then, 
neglecting an insignificant constant, one obtains the lattice potential 

VLix) = Vosm\Kx) (4) 

where k = k cos 9 — tt/A with k ~ 27r/A, A is the laser wavelength, Vq = h^lj^J {ujl — 
f2o) is the amplitude of the OL, = |d • E|//i is the Rabi frequency, and HHq is 
the energy difference between the ground and first excited atomic states. Hereafter 
we neglect losses assuming that the laser beams are far-off-resonant. 



laser beams 




Fig. 1. Typical geometry of a BEC in an OL. Here ki, k2, and the x-axis lie in the same plane. 
|ki| = |k2| = k, |ei| = |e2| = 1 and ei |j 62. 

The BEC dynamics is characterized by one more relevant length scale - the 
healing length ^ = (87rn|as |)~^/^, where the s-wave scattering length and n oc 
N /a\aQ is a mean particle density. Thus we collect four** characteristic scales which 
will define the dynamics of a BEC in an OL: {ao, a^. A, . Alternatively, one can 

''I as I being of order of a few nanometers is much smaller than other characteristic scales. 
''Although, we introduced A through the laser wavelength, bellow we use this notation in a gen- 
eralized sense of a period of the potential. For example, in the case (5) A = 2K(m)/K. 
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consider a hierarchy of the respective energy scales: {huj±,HuQ,Eii,gon}, where 
Er = h^K? /{2M) is the recoil energy. 

It worth pointing out that by means of using several laser beams one can create 
more sophisticated periodic potentials. From the theoretical point of view potentials 
expressed through Jacobi elliptic functions play a special role, because they may 
allow obtaining exact solutions^^ (see also Section 6). In particular, we will be 
interested in the potential 

^hix) = VoST?{k,x I m) (5) 

which is reduced to (4) in the limit m. 0. Hereafter we use the standard^"^ nota- 
tions pq(KX I m), where p=s, c, or d, and q = n, for the Jacobi elliptic functions 
with the parameter < m < 1. Indeed, using the expansion of the elliptic sine^^ in 
terms of the small parameter q = exp ^— tt ^^^^^^ ^ <C 1, where K(m) (and E(m) 
below) is a complete elliptic integral of the first (respectively second) kind, one can 
approximate 



Vl{x) 



2Vo7t\ 



mK2(m)(l -q)2 



TTKX „ 2TrKX 



1 — (1 — 2q) cos — — r- — 2qcos 



K(m) ^ K(to). 

In practice this means that by using four laser beams one can approximate potential 
(5) for a rather large interval of m, ranging from zero up to approximately 0.9 with 
accuracy higher than 99%. 

3. One-dimensional models 

Let us consider a cigar-shaped condensate: a± <^ ao, at a, relatively low density, 
when the energy of two body interactions is much less than the kinetic energy in 
the transverse direction, i.e. when 



e 



2 



Mas 



« 1. (6) 



Then, the system becomes effectively ID^^. In literature there exist several approx- 
imations reducing the original GP equation to an effective ID GP (or NLS) model. 

A self-consistent way, allowing one to consider different relations among parameters 
and to control magnitudes of neglected terms, consists in use of the multiple-scale 
expansion with respect to the small parameter e given by (6). We outline the method 
for a specific relation among the parameters, referring for more details to^^'^^'^'^, 
and list other relevant cases. 

To this end we scale out the variables t = uj±t/2, x — x'/a_L, — r^/a^, 
and the wave function ^ = aj.la^l^/^^', and consider a long trap with a rapidly 
varying OL, what is quantified by the relations a_L ^ A < e^ao. In this case 
GP equation (2) is rewritten as 



= -Ai + fi, £o = -dl + u'^x' + j^Ui^x). 



I 



I 
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where a =sign(as), v — ijJq/uj±^ is the aspect ratio of the trap, U{'yx) = VL{x)/E}i 
is a dimensionless periodic potential varying on a unit scale, U{jx) = U{'y{x + A)), 

A = A/ax^l,and7^ = |ij = 42f. 

Next we consider the linear eigenvalue problems 

^0 ^a,q{x) = £a,q (pa,q{x), Cj_ 0/m(fi) = Elm 4>lm{i±), (8) 

where the indexes a and q stand for a number of a band and for a wave vector inside 
the first Brillouin zone (BZ), and I and m stand for radial and magnetic quantum 
numbers. The operator Co has a discrete spectrum, which however approximates the 
band spectrum in the limit <C e, what justifies the use of the "band" terminology. 
Now (pa,q{x) is approximated by a Bloch function (BF) of the respective periodic 
potential, i.e. by a solution of the eigenvalue problem 

+ U{x)(pa,q{x) = £a,q'Pa,q{x) , (9) 

where taking into account that aj_ ^ A we put 7=1 and provisionally dropped 
tildes. Below we will restrict the consideration to the states bordering an edge 
of the first BZ where q = qo = 27r/A, and respectively will use abbreviated 
notations ipa{x) = (pa,qo{x) and £^"^ = £a,qo- Also we will be interested only 
in evolution of the background transverse state {I = m = 0) corresponding to 
4'oo{^±) = G^p{—f±/2)/^/n and to £^00 = 2.^ All eigenfunctions are considered to 
be normalized to one. 

Now we look for a solution of Eq.(7) in a form 'J' = etjji + e'^tp2 + • • • , where 
ipj {j = 1,2,...) are functions of f_L and of scaled variables tp = eH and Xp = 
e^x, regarded as independent. It follows from (3) that Vi is normalized as follows 
JlV'ipd^f = ^ what means that = Then the leading order of the 

solution can be searched in a form of the modulated ground state 

^1 = ^^(a;i,t2)e-^(^'"'+2)*°e-^-/Va(a;o). (10) 

Here A{xi,t2) (as well as Bn{xi,t2) below) is a slowly varying envelope amplitude, 
where by convention we indicate only the most rapid variables (i.e., for example, 
-4(xi,t2) means that A depends on all {0:1,2:2, • • •} and {^2,^3, ■ • •})• 

After substituting the expansions into (7) one collects the terms at each order 
of e (for the set of so obtained equations see e.g.^^'^''). The first order equation is 
satisfied by ansatz (10). In the second order one obtains 

= ^e-*=i/2e-'(^'"'+2)*<' V B^,^{xuti)<fia'{xo) (H) 

where B^^a = g(oirig(o) ^a^i^ and Ta-,a = ~ Jq ^ai ^^xq. Thus V2 always exists 
when the ground state is excited (see Fig. 2). That is why, following the nonlinear 



■^Strictly speaking intcr-lcvcl transitions are not forbidden. They however can be originated by 
lattice perturbations or by initial conditions which will not be considered here. 
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optics terminology^*, we call ip2 a companion mode. Also it follows from the second 
order equation that dt^A = 0, i.e. A = A{xi,t2)- 

Eliminating secular terms from the third order equation we obtain 

idt,A + i2M„)-^dl^A - Xa(r\A\^A = 0, (12) 

where the effective nonlinearity, Xai and the inverse effective mass, Ma, are given 

by 

Xa = 2f\^a{xordxo, and -L = 2 + 2 J] ^^^^f^ . (13) 

By substitutions A \/2/xa'4', ^2 t, Xi ^ x, and — > M Eq. (12) is reduced 
to a dimensionless NLS equation 

i5tV = -^52^ + 2a|V'|V- (14) 

To conclude wc have to clarify the meaning of the coefficient M~^. This can be 
done^*, by means of the so-called kp- method, known in the solid state physics^^. 
Namely, it turns out that in the limit ao/A — > oo one has^^ — '^^a,q/d^ ■ 

By analogy one can consider other relevant cases as follows (notice that in each 
case the operator Cq is redefined): 

- Short trap and rapidly varying optical lattice: aj_ ^ A ~ ~ eao- Wave 
function satisfies the equation (all tildes are suppressed) 

*5tV = -^S2V' + z^VV + 2a|^|V (15) 

The eigenfunctions are now the exact BFs solving Eq. (9). 

Long trap and smooth optical lattice: ^ eA ^ < e^ao- Wave function is 
governed by the equation (in the final equation all tildes are suppressed) 

idti) = -dlil) + U{x)i) + 2a\i)\'^'il). (16) 

Now the periodic potential has the scale equal to (or larger than) the healing length, 
and thus excitations in a BEC have also scales of order of the lattice period. It is also 
assumed that 7 = £7, where 7 is the aspect ratio and 7 = (2a_L/A)(ao/A/'|as|)^/^ = 
0(1), and that 72;7(7i) = t^tj(xi). 

- Short magnetic trap and smooth optical lattice: a_L eA ~ eao. Wave 
function is described by the equation 

idtip = -dli) + u^x^ip + U{x)ip + 2a-|V'| V- (17) 

Below we use one of these models depending on the effect we are interested in. 
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4. Modulational instability 

As it is clear from the preceding section, the eigenfrinctions of the linear problem (9) 
play an important role in the theory of low density BECs. When inter-atomic inter- 
actions are neglected, they represent renormalizcd wave functions of single atoms, 
whose evolution in time is trivial. Interactions drastically change this situation, lead- 
ing in some cases to instability of the ground state. This is a phenomenon usually 
referred to as modulational instability (MI) and is well studied in hydrodynamics, 
plasma physics, nonlinear optics, etc.^°. 

In order to develop a theory of MI of a EEC embedded in an OL we, following 
Refs.^^'^'', consider first the model (14) and look for a solution in a form of the 
plane wave V-' = (p + ae*'^^^*^'^^' _|_5g-*(S2t-K;:r)-)j,-ip i ^ -where |a|, |6| p. Linearizing 
(14) with respect to a and b one obtains the dispersion relation = Z{Z + Aap"^), 
where Z = K? /2M, from which it follows that Q. has an imaginary component, and 
thus the plane wave is unstable, if 

Z{Z + Aap^) < 0. (18) 

This condition is satisfied for large wavelengths (small wavcvcctors, |/C| < 2p^J\M\) 
of the perturbation if aM < 0. Taking into account that the effective mass has 
different signs at different edges of the forbidden gap (see the example in Fig. 2) 
we conclude that one of the two matter waves bordering the gap is modulationally 
unstable. 

To explain MI also in a different way we consider long wavelength excitations 
and substitute if) = pexp{—iV,t) in (14). This gives a dispersion relation fl = ap^. 
Recalling that ip is an envelope, the total frequency of the wave is given by £'^°'^+e^Q. 
(see (10) where 2 must be omitted in the ID case). Then, if ct > 0, for the upper gap 
edge f corresponding to the positive effective mass, M > 0, the energy shift due 
to the nonlinearity is positive, i.e. the total frequency falls into the allowed band 
(fds Fig. 2). At the edge £^'^\ where M < 0, the energy is shifted upward in the 
gap (fbs Fig. 2), where plane wave cannot exist, what means instability. If cr < 0, 
the same arguments are applied to the £^^^ and S^'^^ exchanged. 

For smaller wavelengths of perturbation MI can be studied numerically^''. In 
the right panels in Fig. 2 we show the results of such study of the scaled out 3D GP 
equation with the potential 

C/(.x) = Acos(2a;). (19) 

The upper and lower panels show respectively stable and unstable evolution of the 
Bloch states bordering the edges of the first lowest gap. The localized pulses emerg- 
ing from the periodic wave arc the bright matter solitons (sec the next section). 

The fact that an OL can result in the instability of an otherwise stable BEG 
has recently been confirmed experimentally^. 

The above analysis of the MI can be expanded to 2D and 3D periodic 
potentials'^. In higher dimensional cases the inverse effective mass is re- 

placed by the tensor of inverse masses M~\ = dfjdjSa,^ where df} = 
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Fig. 2. Left panel: The band gap structure of problem (9) with A = 1. In a BEC with a posi- 
tive scattering length f^^' (respectively Mi < 0) and £^^^ (respectively M2 > 0) correspond to 
unstable and stable Bloch states. If the main mode, Vi (see (10)), is excited against a carrier 
wave background fiix) (the filled in symbol) then the companion modes ip2 (see (11)) are excited 
with the same wavevector but belonging to different zones (open symbols). The bright soliton 
frequency fbs (the lowest dashed line) is shifted toward the gap. In the small amplitude limit the 
shift is given by rj'^ /{2\M\) (see (21)). The upper dashed line shows the energy of the dark soliton, 
shifted toward the allowed band (in the small amplitude limit the shift is given by r]^/\M\ (see 
(22)). Right panels: Numerical solution of 3D GP equation (7) with cigar shaped geometry and 
ID OL in the x-direction as in the left panel. As an initial condition we use either ipi[x), well 
approximated by \/2sin(3::) (the right bottom panel), or (fi2ix), approximated by \/2cos(x) (the 
top panel), both modulated by the Gaussian function /^-vx /2 



r = (ri,r2,r3) = {x,y,z), and /3, 7 = 1,2,3. When the multidimensional lat- 
tice potential is separable, say when VL(r) = J2^=i cos(Kjrj), the tensor M~^^ 
is diagonal, and the analog of model (14) reads^ 

i^ti^ = E ^ai-^'V' + 2a|V| V- (20) 

j 

Now the dynamics is more rich since BFs corresponding to different spatial 
dimensions contribute to the wave dynamics. The parameter Z, introduced in (18) 
has now the form Z = J2j ^ajj^]- Depending on the curvature of the energy 
surface £a,k in the given spatial direction, /?, the respective band edge can be 
called cither stable, if ^^^^§^(J > 0, or unstable, if — ^S^o" < 0. In Fig. 3 wc show 

unstable evolution of the BF in the two spatial dimensions (results from Rcf.'^^). The 
figure clearly demonstrates emergence of a coherent structure consisting of spatially 
localized matter waves. The symmetry of the emerging structure is defined by the 
symmetry of the lattice and by the characteristic wavevector of the most unstable 



2, 2008 8:6 WSPC/INSTRUCTION FILE BK'review 



9 



mode (see^^ for the details). Existence 
about possibility to generate 2D and 3^ 
indeed been observed in numerical studi 




X 



F such structures leads to a supposition 
matter gap solitons. Such solitons have 




X 



Fig. 3. The example (from Ref.^^) of unstable evolution of BEC atomic distribution in a 2D OL 
aecording to equation (20) with Vi = V2 = 1.0, ki = K2 = 2.0, c = 1/2, L = 14n. (a) The 
initial state; '4){x,y,0) = 0.5 sin(x) sin(j/) perturbed by 6V'(a;, = 0.01 sin(0. 125a;) sin(0. 1252/). (b) 
Formation of soliton-like excitations due to MI at t = 85. 



5. Matter solitons and intrinsic localized modes 

MI instability considered above is not the only important consequence of the exis- 
tence of the efFcctivc mass in the theory. Depending on the sign of Ma an OL acts 
as a focusing medium, if Ma < 0, or defocusing medium, if Ma > 0. Taking into 
account that by using a moving lattice one can obtain either of two signs of the 
effective mass in the case when the nonlinearity dominates dispersion, the lattice 
can be viewed as a compressor of the matter waves^^, which was recently achieved 
experimentally^^ . 

When nonlinearity and dispersion are in balance, creation of stable structures, 
termed solitons is possible. In particular, NLS equation (14) possesses^^ either a 
bright soliton solution 

if Ma < 0, or a dark soliton solution 

" ^/iW] '''''' (22) 

if Ma > 0. The parameter rj characterizes widths and amplitudes of the solitons. It 
is natural to refer to these solutions as matter gap solitons. Well known in nonlinear 
optics of photonic crystals^* matter gap solitons have been predicted in Ref.^^. Very 
recently bright solitons in a condensed cloud of ^''Rb atoms have been generated 
experimentally^^ . 
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As one can see from Eqs. (14) and (15) one of the main features introduced by 
an OL with a small lattice constant is a change of the spectrum. This is expressed 

in terms of the effective mass (13) which can change its sign. Thus even in the case 
of repulsive (attractive) atomic interactions in the presence of OL one can excite 
bright (dark) matter solitons, what is impossible in homogeneous condensates. As 
it follows from the multiple-scale analysis solitons (21), (22) are characterized by 
the following main features: (i) They are low-amplitude, having frequencies close to 
the gap edge (shifted toward forbidden and allowed bands, respectively); (ii) They 
extend over many lattice periods; (iii) They cannot be excited without companion 
modes, what contributes to instability and finite life-time of these entities; (iv) They 
arc not exact solutions of the GP equation, and as such can be deformed during 
propagation what results in their decay with time (for radiative losses of matter 
gap solitons see Ref.^^). 

5.1. Tight-binding approximation 

If a gap, i.e. \£^^^ — £^^^\, is large enough, what happens if the potential rchef is 
deep enough, then when the energy shift \£ — £^"^ | outwards the gap edge £^"^ 
increases, a matter soliton does not disappear but becomes strongly localized on a 
very few lattice cites^^'^^. Now one cannot use model (14) any more, and has to 
explore Eq. (16) describing the situation where the healing length is of order of 
the lattice period. The respective approach is often referred to as a tight-binding 
approximation^^ '^"^ . We introduce it following Ref.^^ where a mathematical back- 
ground for the method was elaborated. 

Periodicity of the lattice potential implies periodicity of the eigenvalues of (9): 
£a,q+qo = ^a,q = ^a,-q whcrc qq = 27r/A what allows One to expand £a,q in the 
Fourier series: £a,q = X^n '^'l"'''*''"^' where LUn,a = ^ Il°q{^/2^a,qe~'"'"^dq. Next, 
we introduce the Wannier functions'^^ (WFs) 



(here n is an integer). Some examples of WFs are shown in Fig. 4a. Taking 
into account that WFs are normalized, J^^Wn^a{3:)wn' ^a'{x) dx = Sn^n'^a,a'j 
make up a complete set, '^na'*^n,a{x)wn,a{x') = S{x — x') and can be cho- 
sen real"^^, one can look for a solution of nonlinear problem (16) in the form 
'ilj{t, x) = Cn^ait)u>n,a{x) where the expansion coefficients solve the set of equa- 
tions 





ai,a2,o:3 711,712,^3 



with the nonlinear coefficients given by 




dx. 



(25) 
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Fig. 4. (a) WFs of the three lowest bands for the potential (19) with A = 5. The WF of the 
first zone (bold solid line; thin solid line and dotted line show WFs of the second and the third 
zones, respectively) coincides in the leading order with the shape of a Wannier soliton. (b) Atomic 
distribution over the sites in a discrete soliton corresponding to the Wannier soliton. 



In its general form, Eq. (24) is equivalent to (16) and thus it is not solvable. 
However, it allows reductions to a simpler lattice in a number of important par- 
ticular cases. Indeed, Wn,a{x) can be chosen to be localized around x = nA and 
exponentially decaying. If the OL depth is large enough, among all the coefficients 
^aalalal those with n = ui = n2 = are dominant and other terms can be 
neglected. Moreover, in this case the lowest zones are very narrow and the Fourier 
coefficients Cjn,a rapidly decay |t!>o,a| ^ ^ |'^n,a|> > !)• Then restrict- 

ing consideration to the only band a one reduces (24) to the discrete nonlinear 
Schrodinger equation^^ (DNLS) 

HC 

i— ^=w(C7„_i+C„+i) + <7|C„|2C„ (26) 



where C„ = c„_q exp(— «a)o,at)-\/2V1^"™" and lu = LUi^a- In the BEC theory DNLS 
(26) is referred to as the tight-binding approximation^^'^"^'^^ . 

Now one can obtain the simplest stable solution which represents strong density 
localization near one of the potential minima, say near n = 0. To this end we restrict 
the consideration to cr > and notice that considering the first zone a = 1 already 
for Vq/Er = -5 we have^^ a)o,i « -2.1152, a>i,i « -0.0192 and u)2,i « 0.0002, and 
thus Un,i can be regarded as a small parameter already for n> 1. Next we fix the 
maximal density |CoP = 1 and look for a symmetric solution, C-n{t) = Cn{t) in the 
form of the expansion: C„ = e~*^^* Xijln where f2 = Y^^^qO.'^^^ Cj^ . Direct 

substitution of this expansion into (26) allows us to compute^'^: O = l + 2uJ^ + 0{Cj^) 
and Cj = e~'^*c2'^ +o(a;^"'""^). Thus the decay of the atomic population is exponential: 
|C„| cx e^''l"l, where 77 = — \tlCj (a thin line in Fig. 4b). For the above values of 
the relation between the potential energy amplitude and the recoil energy about 
Ri 90% of the condensed atoms are concentrated in the vicinity of the central site. 
Taking into account that the shape of a discrete soliton is well approximated by a 
WF, such discrete solitons can be called Wannier solitons^^. 
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5.2. Matter solitons. General approach 

The tight-binding approximation is vahd when the lowest allowed zones are very 
narrow. Then one can approximate a)o,i ~ £^^^ and a)o,2 ~ £^'^^ ■ Hence shift of the 

energy of a discrete (Wannier) soliton toward the forbidden zone is not small. Then 
a natural question about transformation between these gap solitons and discrete 
solitons, when the energy changes, arises. 

To answer this question we following Ref.""^ look for a solution of (16) in the 
form tp{x,t) = e~^^*(fi{x) which leads us to the equation 

fxx + {£- U{x) - 2a\(p\'')>fi = (27) 

where £ plays now a role of a spectral parameter. Next, we separate the amplitude 
u{x) and the phase 6{x) ofip{x): (p{x) = u{x)e~^^'^^\ By substituting this expression 
into Eq. (27) we obtain the link u'^dx = v and the equation for u{x): 

Uxx-:^ + {£- U{x) - 2au^)u = 0. (28) 

The constant v can be interpreted as a normalized component of the condensate 
velocity, because v oc {tp^px — i^i'x)- 

Due to the second term in Eq. (28), u{x) can vanish only for the zero velocity 
V = 0. Thus any stationary solution acquiring zero value at any point, including 
bright matter solitons (they have zero boundary conditions) can be excited only 
with a phase independent on x. Without restriction of generality in what follows 
we consider 6{x) = and thus ip{x) = u{x), where u{x) solves the equation 

Uxx + {£- U{x) - 2au^)u = 0. (29) 

Eq. (29) can be interpreted as a Hamiltonian system with one and a half degrees 
of freedom. Let us, following**^, consider periodic potential (19). Since Eq. (29) is 

invariant under shifts of the period of the potential, one can describe the solution by 
studying the area preserving Poincare map T : {ux{x), u{x)) — > {ux{x + tt), u{x + 
tt)) in the plane {ux,u). The riTr-periodic solutions of Eq.(29) correspond to the 
orbits of period n of the map T. i.e. to fixed points of the map T". These fixed 
points can be either hyperbolic or elliptic^^. Solitons correspond to homoclinic 
and heteroclinic orbits of the Poincare section, i.e. to sequences of points {i?„ = 
«,«„), n e Z,u„ = u{xq + Trn),u'„ = WxU=j,o+^„}, = ^^^^ ^ 

tends to the same (homoclinic orbit) or different (heteroclinic orbit) fixed points 
of T as n = ±oo. Each hyperbolic fixed point O has stable, Ws{0), and unstable, 
Wu{0), manifolds, whose intersection means the existence of a homoclinic orbit. 
If the map T is defined in the entire plane and is smooth, a single transversal 
intersection of these manifolds implies the existence of infinitely many homoclinic 
orbits of T that may be coded by Bernoulli sequences*^. Properties of T depend 
on the sign of a. For cr < 0, the map is well defined in the whole plane, so that one 
expects the existence of infinitely many stationary soliton-like solutions of Eq. (29) . 
On the contrary, for a > 0, the solutions of Eq. (29) may blow up in modulus due 
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to the sign and the power of the nonUnearity, so that the map T cannot be defined 
in the entire plane {ux,u). In this case a number of solitons is drastically reduced. 

Bright matter solitons correspond to solutions of Eq.(29) that vanish as x 
±00. In terms of the dynamical system generated by the map T, they are linked 
to homoclinic orbits of the fixed point O = (0,0), which must be hyperbolic. This 
happens when £ lies in a forbidden zone of (9) . Here we restrict consideration to the 
first forbidden zone: £^^^ < £ < £^'^\ Eq.(9) has only one solution (up to constant 
factors) which decays as a; ^ oo: </?o(a;) = P{x)e^'', where P{x) = P{x + 2^) and 
H the Floquet exponent. In this case, the eigenvalues of the map T at the point 
(0,0) are Ai^2 = e^'^'^ ■ From Eq. (9) it follows that the asymptotic behavior of the 
bright soliton solution is given by the expression u{x) ~ G(fio{x), where G is a real 
constant, what allows numerical constructing the solutions of (29). It was found in 
Ref.**^ that there exist two (up to sign) simple states (even and odd) centered in a 
maximum (minimum) or in a zero of the lattice potential. All other bright soliton 
solutions can be interpreted as bound states of these elementary entities. 

Fig. 5a represents a stochastic layer obtained by iterating the Poincare map T 
starting with an initial point close to (0, 0) (about half a million points). The opened 
circles on top of the stochastic layer arc the points {m^, m„} that correspond to the 
bright soliton solution of Eq. (29), explicitly shown in Fig. 5c. When the energy £ 
is lowered to the bottom of the gap, f^^-', the stochastic layer becomes thin and 
the amplitude of the corresponding solution tends to zero, what corresponds to the 
solution (21) obtained in the weakly nonlinear limit (Section 3). When the energy 
approaches the upper boundary of the gap, £^'^\ the soliton amplitude growth. 
Since the gap width is finite there exist upper bounds for the amplitude of a matter 
soliton and for a number of particles constituting it. 

Matter solitons on a periodic background also can be constructed using the 
dynamical system generated by the map T. They correspond to homo- and hete- 
roclinic orbits of hyperbolic fixed points of the maps T" (n = 1, 2, . . .), other than 
the point (0,0). Consequently, the background must be only a hyperbolic periodic 
solution of Eq. (29). It has been found in Ref.*° that there exists a set of tt- and 
27r-periodic solutions of Eq. (29). In the case 6{x) = 0, one can find solutions that 
are characterized by Fourier expansions with either cosine or sine terms (up to 
7r/2-shifts). 

Each periodic solution of the hyperbolic type can be related to four simplest 

soliton solutions against its background. Two of them are bright and dark symmetric 
solutions, i.e. having a maximum or a minimum of amplitude at the center (the 
former are also referred to as anti-dark solitons), and the two others are bright 
and dark antisymmetric solitons. In the case cr < 0, we have observed that all the 
four soliton solutions may coexist for some range of the values of £, A, and the 
amplitude of the background periodic solution. They may also bind together and 
form complex multisoliton states. The situation is different for tj > 0. In this case, 
there exist two kink-shaped solitons, which can be treated as transitions between 
two periodic backgrounds of different sign. An example of numerical study of the 
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Fig. 5. Phase portraits obtained by iterating T for Eq. (29) with the parameters: (a) A = — 1, 
£ = 0.55, (7 = 1/2 and (b) A = 1, £ = 0.3, a = —1/2. (c) An example of a bright soliton 
corresponding to open circles in (a); Examples of a bright (d) and dark (e) solitons against a 
periodic background (squares and circles in (b) correspondingly). 



Poincare map T corresponding to existence of a dark and anti-dark solitons is 
shown in Fig. 5b, d, and e. These solutions are counterpart of bright and dark 
solitons against a background in the theory of DNLS equation^^ (26). 

In the present text wc do not discuss stability of matter solutions. In the case 
of small-amplitude solitons the evolution equation is obtained by excluding secular 
terms and thus guarantees stability in some sense. In a more general case the 
stability can bo established by the direct numerical simulations or by providing the 
linear stability analysis^"^'^^. Wc also notice that the dynamical system approach can 
be applied directly to a lattice approximation to the GP equation^'"'. Meantime, a 
complete classification of all available matter soliton solutions like this was recently 
donc'*^ for the discrete intrinsic localized modes of the DNLS model remains to be 
an opened problem. 
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6. Nonlinear periodic waves 

As it is well known, NLS equation (14) possesses exact periodic wave solutions 
which, also like solitons, can be obtained by means of the inverse scattering method. 

In our case this means that one can excite not only localized matter waves hut 
also nonlinear periodic matter waves. Periodic waves survive in the limit when 
the healing length becomes of order of the lattice constant, i.e. in the case of the 
evolution equation (15), what has been shown in Section 5. A remarkable fact is 
that for some particular potentials one can find exact analytical periodic solutions 
valid for the whole range of the relation between the lattice constant and the healing 
length. This is the case of potential (5), studied in detail in^^ and briefly summarized 
in the present section. 

Let us focus on stationary solutions of Eq. (16). In contrast to bright matter 
solitons, periodic waves may do not acquire zero values on the real axis, and then 
they may posses .;:-(iependcnt phases (the so-called nontrivial phase solutions). If 
the phase is constant (then it can be put zero) the solution is referred to as a trivial 
phase solution, and it solves Eq. (29). Here we discuss only the last case. 

We observe that if stationary periodic solutions of the NLS equation (14) with 
M = 1/2 arc known, wc denote them -0 = e~*^*wo(a;), one can easily verify, that 
the function ip = e^'^*?i(.T) with u{x) = Y^Mo(a;), x > Oj solves (16) with a 
periodic potential U{x) = 2(t{1 — x)\uoix)\^. Alternatively, Eq. (29) allows one 
to construct a lattice potential U{x), in an explicit form, when the desired wave 
field Uo{x) is given, subject to the requirement of a regular behavior of uq^xx/uq: 
U(^x) = "^^^ + £ — 2c7Uq. In particular, one can verify that this is the case of the 
Jacobi elliptic functions pq(«;a;,m), where p=c,s or d, and q = n (more generally 
one can consider, say, Upq(x) = (j4pq'^(K.x, m ) + S)" with B >\A\,a and /? being 
real integers). Respective solutions, which subject to proper choice of £ lead to 
potential (5), are as foUows^^: 



ij,n{x,t) = e''^'(i+'")*u,„(a;), = ^ (^n^m - asn{KX \ m), (30) 

Vcn(a;,i) = e^['^'(2r„-i)-^''l*wen(a;), u^x) = ^(y-K2m^(7cn(Ka; | to),(31) 

i^dn{x,t) = e'^'^'^^-^^-^Kd^x), Ud^{x) = ^j~(^^^adn{Kx\m), (32) 

where the parameters k, k > 0, and m, < m < 1, must be chosen to make an 
expression under the square root positive. 

In the absence of the lattice potential, Vq = 0, the sn-wave exists in a BEC 
with a positive scattering length, a — 1, while cn- and dn- waves exist in BECs 
with negative scattering length, a = —1. They are solutions of the NLS equation 
(14) with M = 1/2. The trap potential changes the situation, allowing existence of 
sn-waves (cn- and dn-waves) in BECs with negative (positive) scattering lengths. 
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If one considers the limit rn —> 1 at Vq = (or other allowed values of the 
potential amplitude), one obtains the transitions from the periodic waves to the 
solitons, given by (21) and (22) with Ad = 1/2 and k = -q: ipsn V-'ds, "^cn ^ i'hs, 
and -i/'dn — * V'bs- This means that by manipulating the parameter m one can provide 
transformation among different solutions (see Section 7). 

Finally it worth to mention that the linear case, corresponding to very low 
densities, i.e. to the limit of zero amplitudes, is achieved when m what is 

possible when Vq > Q. In this case, the above solutions are transformed into the 
particular Bloch functions of the potentials Vl = Vosv? (^kx \ 

7. Management of matter waves by Feshbach resonance 

When a periodic wave is given by one of the Jacobi elliptic functions one can provide 
deformation of this wave to a sequence of localized matter waves by modulating pa- 
rameters m and k. On the other hand, in experimental*'^^ and theoretical*^'*^ works 
it was shown that an effective tool for manipulating BECs is the Feshbach resonance 
(FR), allowing one to change the scattering length by varying external magnetic 

field B{t) near the resonant point Bq (Ref.^''): a.,(f) = as{0) ^1 + ■^^^(tj^ (here 
as{0) is the asymptotic value of the scattering length far from resonance, and A 
is the width of the resonance). Since change of the nonlinearity results in change 
of m (see below), trains of solitons can be generated from a periodic wave in a 
controllable manner^". 

Let us consider the problem of soliton generation in more details (generalizing 
earlier results^" to the case of a BEG in an OL). Temporal dependence of means 
the dependence on time of the nonlinear coefficient: go — > goi^'^t) = ga ■ .9(^2) 
where we assume slow dependence on time and normalization g[0) = 1. Then one 
can repeat the arguments of Section 3, substituting go by 50(^2)- The modification 
of each of the evolution equations (14) (17) is reduced to the modulation of the 
nonlinearity by the function g{t). In particular, (16) is reduced to 

+ 'ipxx - U{x)ij - 2aff(t)|V'| V = 0. (33) 
g{t) is a positive function (the only case considered below) by the ansatz 

V(a.,i) = ^%^ = e-fW/M^,C), (34) 

where C = C(i) = 2 J*j{t')dt' = In g{t), C(0) = 0, and v{x,0) = Jp{x,Q) Eq. (33) is 
transformed into a dissipative NLS equation 

ivt + Vxx — U{x)v — 2a\v\'^v = ij{t)v (35) 

where 7(f) — gt/2g = Ct/2- Thus, for slowly varying g{t), \^{t)\ <C 1 and the right 
hand side of Eq. (35) can be considered as a small perturbation. 

Let us focus on perturbed evolution of periodic waves assuming that initially 
they are given by v{x,0) = Upq{x) where Upq{x) is one of the functions given by 
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(30), (31) or (32). Under influence of the dissipative perturbation a wave shape is 
changed what can be described by variations of the parameters m and k assuming 
them to be slow functions of time t (the adiabatic approximation). To this end wc 
observe that the coefficients in the dissipative model (35) are independent of a;, and 
thus the perturbation does not result in change of the period A of the nonlinear 
wave what can be expressed as follows 

-^A(m, k) = 0. (36) 

Next it is a straightforward algebra to ensure that the integral N = \v\'^dx is 
governed by the equation 

dN(m,K) , ,„„, 

— ^^=iV(m,«). (37) 

If the expressions for A and N in terms of m, and k arc known, Eqs. (36) and 
(37) reduce to a system of two differential equations of the first order for m and k. 
The form of this system depends, of course, on the choice of the initial pq-wave. 
Below, we concentrate on two examples of transformations of periodic waves in 
trains of bright and dark solitons (for more details see Refs.^^'^°). 

7.1. dn-wave in a BEC with a negative scattering length. 

Consider a BEC with a negative scattering length, a = —1, which has initial form 
of the dn-wave u^nix), given by Eq. (32). We compute A = 2K{m)/K, and N = 
{2K'^m-Vo)E{m)/{Km), and from Eqs. (36), (37): 

dn K(2K2m-Vo)^. .r-c^^ . 

— = ^ ^E(m) [E(m) - m'K(m)] , 

ac, A /Qo\ 

dm 2{2h-n} - \\{)inm' 

— = E(m)K(m) 

where A = Vo [E'^{m) + m'K^{m)] + 2k^m [E'^{m) - m'K^{m)] and m' = 1 - m is 
a standard notation. 

In Fig. 6a we present a numerical solution of the system (38) which illustrates 
creation of a train of bright solitons in the minima of the OL. 

To check the adiabatic approximation we solved numerically (Fig. 6 b, c) the di- 
mensionless 3D CP equation for cigar shaped condensate with cylindrical symmetry 
and a = 1: 



idt^ = -{dl + ^drrd^ V + K.xt(x,r)V' + 2e*/nV'lV 



(39) 



where = + and Vexi{x.,r) = + //^x^ + VbK^sn^(Ka; | m). We switch off 
the FR at Tq and in the further calculation we take g = <?(To). One again can 
clearly see that after some time a set of bright pulses localized in the minima of 
the OL potential arises. These pulses represent solitons, which would survive for a 
sufficiently long time after the lattice potential is switched off. 
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Fig. 6. (a) Adiabatic modulation of a dn-wavc solution with m(0) = 0.11 and k(0) = 0.75 at = 
(thick solid lines) and at = 3 (thin solid lines) in potential (5) (dashed lines) with Vo = 0.1. (b) 
Initial (t = 0) and (c) final (t = 6 ms) density profiles obtained by numerical solution of (39) with 
fc" = 0.02 and N 10'* atoms of ^hi. Initially the dn-wave has the same parameters as in (a) but 
now modulated by the Gaussian function exp(ya;^/2 + r-^/2). The amplitude of OL is Vb = 0.2^?^. 
The size of the image is 200 x 80 ixm. The chajracteristic times of FR are t = 2 ms and To = 3 ms 
(see the text). 



Similar results of generation of trains of bright solitons in a BEC with repulsive 
interactions can be obtain by considering cn-wave (31) as an initial profile of the 
condensate. 

7.2. Evolution of a sn-wave in a BEC with a positive scattering 
length. 

As a second example we consider deformation of the sn-wave (30) in a BEC with 
a positive scattering length, cr = 1 affected by the FR, which leads to creation of a 
train of dark solitons. Following the steps of the adiabatic approximation now one 
can compute A = 4K(to)/k, N = 2{2K^m — Vq) [K(to) — E(m)] /{Km) and deduce 
from (36), (37): 

dn K{2K?'m — Vq) 



dC A 



[E(m) - K(to)] [E(to) - m'K(m)] , 



dm 2mm'{2n^m-VQ) 

— — -K(m) [E(to) - K(m)]. 



(40) 



dQ A 

Here A = Vb \E'^{m) - m'K^{m)\ + 2^2^ |[E(m) - K{m)f - TOK2(m)|. 

The left panels of Fig. 7 show BEC dynamics in the adiabatic approximation, 
where emergence of dark holes, which can be associated with dark solitons can be 
observed. Notice that in the cases (a) and (c) dark solitons are created in the minima 
and maxima of the lattice potential, respectively. A train of darks soliton can also 
be obtained in the 3D simulations of the model (39) starting from the sn-wave 
as initial profile modulated by Gaussian fmiction cxp(i^.t2/2 + r"^ /2). Like in the 
previous case varying nonlinearity results in creation of solitary pulses. However, in 
contrast to the case of bright solitons, if the parabolic potential is switched off, the 
condensate constituting the background for dark solitons, start to expand rapidly 
passing to the linear regime. 
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Fig. 7. Left panels: Adiabatic dynamics of sn-wave solutions at C = (solid thick lines) and 
f = C/in (solid thin lines) in the periodic potential (dashed lines) for (a) Vo = 0.05, and C/in = 5; 
(b) Vq = 0, C/in = 6, and (c) VJ) = —0.05, C/in ~ 6- Right panels: (a) Initial (t = 0) and (b) final 
(t = 0.8 ms) density profiles of the condensate embedded in ID OL (5) with strength Vo ~ —0.2Eji, 
aspect ratio is u = 0.05 and characteristic time of FR is r 0.06 ms and Tq ~ 0.45 ms. Initial 
profile contains N ^ 5x lO** ^'^Rb atoms. The size of the images is 150 X 40 /xm. (c) Crossection of 
the density |i/!(r = 0,a;)p at initial (thin line) and final (thick line) times corresponding (a) and 
(b). In all cases m(0) = 0.11 and k(0) = 0.75. 



8. Bloch oscillations 

In the preceding sections it was shown that OLs can be employed to create solitary 
and periodic matter waves in a controllable manner. In this and next sections we 

concentrate on the use of OLs for manipulating matter solitons. We start with 
the phenomenon very well known in the solid sate physics^^ - Bloch oscillations - 
consisting in periodic motion (in direct and reciprocal spaces) of a quantum particle 
in a periodic potential. Since an atom in an OL undergoes Bloch oscillations^"'^, BEG 
must also undergo Bloch oscillations when loaded in an OL, what has indeed been 
observed experimentally^'^'^^ and studied numerically^^ In these studies a BEG 
was considered as an atomic cloud, where each atom undergoes Bloch oscillations. 
However, a priori it is not evident whether this phenomenon can be observed with 
matter solitons (considered as single entities). Indeed, as it follows from the effective 
mass approach (Section 3), when in reciprocal space the condensate passes from one 
BZ edge to another, the effective mass changes its sign, and thus a bright matter 
soliton existing at one edge becomes a dispersive wave packet at another edge''^''''^. 

It turns out, that Bloch oscillations, although in a somehow modified sense, 
can indeed be observed^^. Mathematically this possibility is explained by two facts: 
lattice solitons subject to a linear force undergo Bloch oscillations^^^^^ and strongly 
localized matter solitons in the tight-binding approximation are well described by 
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the DNLS equation (see Section 5). 

To illustrate the phenomenon we assume that a BEC is loaded in a vertical OL, 

like in the experimental setting®, and thus besides the lattice potential, a gravita- 
tional field, is applied to the condensate. Alternatively, one can consider an acceler- 
ating lattice, where the linear a;-dependent force appears due to the acceleration^^. 
The gravitational force results in the term mgx"^ (notice that now x is a vertical 
axis, and g is the gravitational acceleration) in the GP equation (2). If one also con- 
siders relatively large amplitude of the lattice potential preventing Landau-Zener 
tunneling between bands^, then in the tight-binding approximation (26) a new term 
gncn^a, where gn = ^^^^^^Lj±g, is to be added: 

dC 

'~dt' " ^ ^"-^^ ^ ^ (^\Cn?C„. (41) 

The obtained system was studied numerically in Ref.^^, where it has been discovered 
that the bright solitary pulse undergoes periodic motion, shown in Fig 8. 




Analytically, the effect can be understood in approximation as follows. In- 
troduce momentum P = X^„Pn and energy E = ^" soliton, where 
Pn = {CnCn-i - C„_ i C„) /(2?) and e„ = (C„C„_i + C„_iC„)/2. It is a straight- 
forward algebra to ensure that 

^=gE + aY,{\Cn?-\Cn-l?)en, 

dF " (42) 

— = -5P-a^(|C„|2-|C„_i|2)p„. 

n 

Thus if populations of the neighbor sites of the OL, defined by |C„p, arc not very 
dispersed | |C„p— |C„_ip |<C g, one can neglect the last terms in the both equations 



I 



I 



2, 2008 8:6 WSPC/INSTRUCTION FILE BK'review 



21 



(42), and obtain that P and E periodically vary with the frequency g and with the 
"dispersion relation" P'^ + E'^ = Pq + Eq where Pq and Eq are the initial values of 
the momentum and the energy of the soliton. 



9. Matter solitons in inhomogeneous lattices 

So far we dealt with homogeneous lattices. It turns out that by modulating an OL 
one can efficiently manage matter solitons as they accelerate, decelerate, oscillate 

or undergo the reflection depending on the type of the modulation introduced^^. In 
particular lattice modulations can simulate attractive and repulsive defects in the 
gap soliton dynamics. 

Let us consider two examples, of a BEC with a negative scattering length, 
(7 = — 1, where the periodic potential is given by one of the expressions as follows 

Vl{x) =\i± e-^'^'(^-^^)'j cos(2x). (43) 

Fig. 9 shows the both types of the lattice modulations ((a) in left and right 
panels) and the respective soliton dynamics (left panel (c), (d) and right panel (b) 
and (c)). In both cases we start with numerically obtained small-amplitude bright 
solitons which exist in the vicinity of the upper band gap edge S^"^^ with S^"^^ —£ <0 
(see Fig.2 and also Section 5.2) ensuring that the soliton has a small initial velocity. 




Fig. 9. Left panels: (a) Periodic potential with defect given by Eq. (43) with upper sign and 
parameters e = 0.05 and Ax = 100. (b) The upper edge of the lowest gap (solid line) and the 
energy of the soliton with the initial velocity f = 0.1 (dashed line). Intersections of these lines 
represent the reflection points, (c) Initial profile and (d) profile of the reflected soliton at time 
t = 150. Right panels: (a) Periodic potential with the defect given by Eq. (43) with lower sign and 
corresponding upper edge of the first band gap (here the dashed line is initial soliton energy). 
Parameters of the defect are the same as in left panel. Profiles of the soliton with different initial 
velocities v = 0.1 (b), and v = 0.3 (c) at times t = 300 and t = 190, respectively. 



The left set of panels of Fig. 9 shows reflection of a matter soliton by the defect. 

This phenomenon can be interpreted as a Bragg reflection as it is explained in 
Fig. 9b. Indeed, smooth modulation of the lattice results in a deformation of the 
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gap edges in such manner that the upper gap edge has bigger energy {S^"^^ + AS) 
than the energy of the gap edge of the unperturbed lattice, f (see Fig. 9b). Thus 
a sohton moving from the left to the right faces the boundary of the gap where its 
velocity becomes zero. The soliton is reflected from the defect, which can be viewed 
as a repulsive impurity. 

In the case of a local decrease of the potential depth, the forbidden gap is 
locally shrinks, the situation is described by the lower sign in (43) and illustrated 
in the right panels of Fig. 9. In this case, in the region of the defect the soliton 
energy is effectively shifted from the gap edge S^'^^ toward the allowed band (see 
Fig. 9a), where the group velocity of the carrier wave has local maximum. Thus the 
modulation described by the Eq. (43) with lower sign acts as an attractive defect, 
and some number of atoms could be captured by such a defect, what happens when 
an initial kinetic energy of the condensate is small enough. This is exactly what we 
observe in numerical simulations shown in Fig. 9b (right panel) . The higher velocity 
matter waves pass through the defect without substantial changes (see Fig. 9c). 

To conclude this section we notice that the above examples (for more details and 
other examples see Ref.^*) show that matter solitons in OLs can be manipulated 
by affecting the carrier wave background, something not possible to do with matter 
solitons in homogeneous BECs. 

10. Conclusion 

In this brief review we have shown a diversity of dynamical effects which can be 
observed in BECs loaded in far-off-resonant OLs, the latter representing an efficient 
tool for managing matter waves. In particular, OLs allow one to obtain solitary 
and periodic matter waves, which cannot exist in homogeneous BECs, to deform 
solutions, to arrange their motion in an a priori given manner, and to manipulate 
matter soliton dynamics by affecting the dynamics of the carrier matter waves. 

Limited volume of the present work however, did not allow us to cover a num- 
ber of other important issues related to BECs in OLs, even limiting considera- 
tion to the mean-filed approximation. Among them we mention matter waves in 
multicomponent^^ and in atomic-molecular^" condensates, the effect of Landau- 
Zener tunneling^'^^ , superfluidity^^, delocalization of matter waves^^, dispersion 
management of matter solitons in OLs'"'"'''®^, generation of matter shock- waves in 
OLs^^, propagation matter waves through finite lattices^^, etc. Finally, we have not 
discussed the relevant stability properties of periodic and localized matter waves. 
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